NACA RM AW321 


f m - A&€Bai 


CA3E FILE! 

COPY 

NACA 

RESEARCH MEMORANDUM 


LIGHT DIFFUSION THROUGH HIGH-SPEED TURBULENT 

BOUNDARY LAYERS 

By Howard A. Stine and Warren Winovich 

Ames Aeronautical Laboratory 
Moffett Field, Calif. 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

WASHINGTON 

May 25, 1956 

Declassified February 26. 1958 


WACA RM A56B21 


NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

LIGHT DIFFUSION THROUGH HIGH-SPEED TURBULENT 
BOUNDARY LAYERS 

By Howard A. Stine and Warren Winovich 
SUMMARY 


The optical transmission characteristics of turbulent boundary layers 
in air on a flat plate with negligible heat transfer were measured photo- 
metrically for ranges of Mach number from 0.4 to 2.5. Free-stream densi- 
ties and boundary- layer thicknesses ranged, respectively, from 0.12 to 
0.93 standard sea-level atmospheres and from I-I/2 to 3-l/2 inches. 

It was found empirically that the scattering from a collimated beam 
of white light which penetrates a turbulent boundary layer depends mainly 
upon the integral across the layer of the difference between the free- 
stream density and the local boundary- layer density. The radiant power 
scattered thus appears to be dependent both upon density fluctuations 
proportional to the change in mean density across the bo\indary layer and 
to an integral scale of the density fluctuations proportional to the 
boundary- layer thickness. The scattered light is deviated in all cases 
through very small angles, the maximum in these tests being about O.OOO6 
radian measured from the direction of primary propagation. The distri- 
bution of energy in the scattered field depends mainly upon the integral 
scale of the density fluctuations and is in excellent agreement with a 
theoretical prediction based on the scattering cross section of Booker 
and Gordon. Comparisons of scattered patterns through two boundary layers 
with those through four show small deviations from Lambert’s exponential 
law. These differences are attributed both to secondary scattering and 
to diffuse refraction at the boimdary- layer free-stream interfaces. 

The results show that significant deterioration in resolving power 
can be sustained by optical imaging devices which receive radiant energy 
through compressible turbulent boundary layers. For example, the included 
angle between two equally bright point objects which can just be resolved 
by a 2-l/2-inch objective is about 2 seconds of arc. If this objective 
looks through a 1-3/4- inch-thick tiirbulent boundary layer at a Mach number 
of 2.5 and a density altitude of 45,000 feet, this angle increases to 
about 8 seconds of arc. 
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The results, moreover, show that photometric measurements in the 
radiation field produced hy the interaction of a plane light wave with 
a tirrhulent boundary layer can be used in conjunction with scattering 
theory to deduce average values of the integral scale and the intensity 
of turbulent density fluctuations. 


INTRODUCTION 


Due to the process of scattering (or diffusion), the quality of 
optical images received aboard an airborne vehicle traveling at high speed 
can be expected to deteriorate if the radiation traverses a turbulent 
boundary layer. The subjective effect of this deterioration can be likened 
to the common sensation brought about when one views a distant object over 
a hot chimney and observes a general bl\irring of outline and disappearance 
of fine detail. It is obvious that airborne optical devices, such as 
reconnaissance cameras, star-tracking telescopes, and fire-control sights 
can be similarly affected by high-speed turbulent boundary layers, and 
thus will suffer a loss in resolving power. 

A qualitative insight into the optical scattering ability of tiirbulent 
boimdary layers can be gained by reference to figure 1 which shows two 
shadow photographs of identical bodies in the Ames supersonic free flight 
wind tunnel (ref. l) at comparable Mach numbers and Reynolds numbers. 

Both photographs were obtained with an electric spark illumination of 0.5 
microsecond, and differ in that the body of figure l(a) was laimched 
through still air whereas the body shown in figure l(b) was launched into 
a counterflow at a Mach number of 2. Thus, the light which exposed the 
film of figure l(b) penetrated two turbulent boundary layers on the wind- 
tunnel side walls, and therein was refracted by the turbulent density 
fluctuations in such a manner as to cause the backgroixnd to take on the 
mottled, grain-like pattern shown. Because of the great difference between 
the speed of light and the speed of the turbulent eddies and because the 
spark duration was sufficiently short , the turbulent motion was "stopped” 
in the pattern shown. However, dinring a finite time interval the pattern 
will, of course, change and each ray of light which enters the tiirbulence 
at a given point will illuminate, in the covirse of time, a finite area 
upon emergence. Thus, a one-to-one correspondence between entering and 
emerging light rays is lost, because, in effect, the turbulence introduces 
a random noise field in the primary beam. 

Using linearized theory of geometrical optics, Liepmann (ref. 2) has 
estimated the rms diffusion angle of a ray for various flight conditions. 

A more detailed scattering analysis based on the electromagnetic theory 
of radiation and the theory of isotropic turbulence has been advanced by 
Booker and Gordon, and refined by Villars and Weisskopf (refs. 3 and k) . 
Baskins and Hamilton (refs. 5 and 6) have measiared the transmission 
characteristics of turbulent boundary layers for a limited range of flow 
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conditions. The results of these theoretical and experimental studies 
show that a significant loss in resolving power can be sustained by opti- 
cal imaging equipment receiving radiant energy which has penetrated a 
high-speed turbulent boiuidary layer. 

The purpose of the present experimental investigation is to provide 
additional data on the optical transmission characteristics of turbulent 
boundary layers over a range of Mach numbers from 0.4 to 2.5, free-stream 
densities from 0.12 to 0.93 standard sea- level atmospheres, and bo\mdary- 
layer thicknesses from l-l/2 to 4-l/2 inches and to correlate the observed 
optical characteristics with the appropriate boundary- layer parameters. 


NOTATION 


C Gladstone-Dale constant (O.llTO ± 0.0002 cu ft/slug for air and 

5200 A light) 


d 


focal plane stop diameter, ft 


D entrance pupil diameter, ft 

E radiant intensity averaged with respect to time, watts/sq ft 

F 

f/no. focal ratio, — , dimensionless 


F 


focal length, ft 


I time-averaged radiant intensity at the focal plane of telepho- 

tometer, watts/sq ft 



I 


integral scale of density fluctuations, ft 


M Mach niamber, dimensionless 


N index of refraction, 1 + Cp, dimensionless 


n number of boundary layers traversed 

P radiant flux (power), watts 

p pressure, pRT, Ib/sq ft 
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gas constant, 1715 sq ft/sec^ for air 

radius normal to optical axis of telephotometer with origin at 
the principal focus, ft 

Tt n 

temperature, , R 

1 + (7 - 1)mV2 

air velocity, ft/sec 

distance through boundary layer normal to surface, ft 
attenuation coefficient, l/ft 

boundary- layer correlating parameter (eq. (25)), dimensionless 
specific heat ratio, 1.4 for air, dimensionless 
boundary- layer thickness, ft 

specific inductive capacity, dielectric constant, N^, for air, 
dimensionless 

scattering angle measured from direction of primary propagation, 
radians 

angular aperture of telephotometer, tan“^ — , radians 

2F 

radiation wave length, ft 

effective radiation wave length, 1.705X10’^ ft 
dimensionless radius normal to optical axis of telephotometer, 

FA 

air density, slugs/cu ft 

density of NACA standard sea-level air, O.OO2378 slug/cu ft 

scattering cross section per unit volume, sq ft/cu ft steradian 

angle of incidence measured from normal to flow direction, deg 

angle between direction of electric vector and scattering 
direction 0, radians 


solid angle, steradians 
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Subscripts 

A absorbed 

B boundary layer 

F flow 

M measured 

WF no flow 

n nijmber of boimdary layers, for example, 1, 2, or 4 

o incident 

P polarized 

R Rayleigh 

S scattered 

sep separated region of boundary layer 

t stagnation 

U unpolarized 

00 free stream 


ANALYSIS 

Attenuation of Light in Turbid Media 


Aggregations of particulate matter, such as air, smoke in air, liquid 
emulsions, colloidal suspensions, etc., are known to attenuate a beam of 
parallel monochromatic light in close accord with Lambert's exponential 
law for homogeneous media (refs. 7, 8, and 9). This law states that 
layers of equal thickness attenuate equal fractions of the incident inten- 
sity. In differential form the law is written 

^ = -a dy ( 1 ) 

The attenuation per foot of penetration, a, is a function of the vibration 
frequency and state of polarization of the incident light, the structure. 
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both atomic and macroscopic, of the particles, and, above all, of the 
degree of chaos (random fluctuations) in the particle nmber density. 
Ordinarily, a is considered to be made up of two parts, an absorption 
coefficient, a^, and a scattering coefficient, ag, to distinguish between 
the Intensity decrease due to disappearance of visible light into radia- 
tion at other frequencies, such as heat, and that due to deflection of 
visible energy from the direction of primary propagation. It is the 
latter which is affected by turbulent density fluctuations in air . 

In the application of equation (l) to the boundary layer, the attenua- 
tion per foot due to scattering, ag, cannot be taken as a constant, but 
must be considered as a function of position because the boundary layer 
is an inhomogeneous medium. However, insofar as a turbulent boundary layer 
can be considered an assembly of particles, it might be expected that 
equation (l) would describe the optical transmission characteristics of 
a series of identical boundary layers, each having an average scattering 
coefficient o,g and thickness &. Such a definition is: 

= ( 2 ) 


where n denotes the nimiber of boundary layers traversed. Equation (2) 
forms a starting point for the experimental procedure used in the present 
investigation because it provides a method for deducing, from the results 
of observations through a multiplicity of identical boundary layers, the 
transmission characteristics of a single boundary layer. In particular, 
the relations between four, two, and one identical boundary layers are 



It should be noted that this result is due to the fxmctional form of 
Lambert’s law and that nothing need be known about the scattering coef- 
ficient and the thickness of the scattering medium, so long as the several 
media are identical and contiguous. 



Attenuation of Light by Isotropic Density Fluctuations 


Lack of knowledge of the turbulence structure in a compressible 
turbulent boundary layer precludes the calculation of ag in general. 
However, the attenuation of electromagnetic waves due to scattering by 
turbulence can be calculated for a tiirbulence model characterized by 
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isotropic density fluctuations having an exponentially decaying auto- 
correlation fimction. The calculation utilizes the scattering cross 
section for this model deduced hy Booker and Gordon (ref. 3)^ 'tke turbu- 
lence structure being defined solely by two parameters^ namely^ intensity 
and integral scale of the density fluctuations. 

One might argue against the worth of such a calculation on two 
counts. First, the turbulence model of Booker and Gordon is physically 
unrealistic because the exponentially decaying autocorrelation function 
implies that the turbulence has no microscale. Second, a two-parameter 
model of isotropic turbulence cannot be expected adequately to describe 
the anisotropic inhomogeneous turbulence anticipated in compressible 
turbulent boundary layers. The extent to which scatter propagation is 
influenced by these discrepancies between the turbulence model of Booker 
and Gordon and more refined models involving three or more parameters 
has been discussed by Wheelon, Muchmore and Wheelon, and Staras in 
references 10, 11, and 12. Whereas these investigations show that the 
Booker and Gordon scattering cross section is, of course, altered by 
choice of correlation fimction and considerations of anisotropy, the abso- 
lute changes in cross section are small for large arbitrary changes in the 
turbulence model. One is therefore encouraged to believe that the rela- 
tively simple calculations based on the results of the Booker and Gordon 
analysis can, at the very least, point out the significant variables and 
aid in the analysis of data. 

The radiant flux (power) scattered by isotropic density fluctuations 
per unit solid angle, per unit incident intensity, per unit macroscopic 
element of volume as given by Booker and Gordon is 


This equation states that the fraction of incident radiant power scattered 
into unit solid angle per foot of penetration through the turbulence, Up, 

is directly proportional to the mean square fractional fluctuation of the 
specific inductive capacity of the medium e (and thus to the density 
fluctuations);^ directly proportional to the cube of the characteristic 

^Using definitions given in the NOTATION, which apply to dry air: 


Up - 



sin2 X 


(4) 





for Cp < < 1 


2 
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size of the density "lumps" (integral scale of turbulence), 1; and 
inversely proportional to the fourth power of the radiation wave length, 
A. The dependence of the scattered power on angle is contained both in 
the denominator of equation (A), where the polar angle, 0, is measured 
from the direction of primary propagation, and in the term sln^’X in 
the numerator, where X is the angle between the direction of the polar- 
ized incident electric vector and the scattering direction 0 . Equa- 
tion (4) can be modified so as to be applicable to natural (unpolarized) 
incident light by replacing the term sin^X by the term l/2(l + cos^q) 
(ref. 13 ) . Then 



(2rtZ)^ 



1 + cos^e 


(5) 


From equation ( 5 ) it can be seen that, in addition to being symmetrical 
about the propagation axis, the power distribution in the scattered field 
for an mpolarized plane incident wave is strongly dependent on the com- 
parative magnitudes of the integral scale I and the wave length X. 

In fact, when I is very much greater than X, which appears to be the 
relevant case for visible light and turbulent bo\mdary layers, and since 
the half -power point of the scattered beam makes an angle of X/4«Z 
with the propagation axis, the scattering takes place mainly in the 
forward direction. 

Using equation ( 5 ), one further concludes that a layer consisting 
of isotropic density fluctuations of macroscopic thickness dy and unit 
cross-section area scatters the fraction of incident power 


agdy = 



ay(w,y)dto 


<iy 


( 6 ) 


Thus, for negligible absorption, the attenuation suffered by the primary 
beam due to scattering in penetrating the isotropic turbulent medium a 
distance y becomes (eqs. (l) and ( 6 )): 

^The analysis of Villars and Weisskopf (ref. 4), which employs a 
turbulence model based on Kolmogoroff • s concept of local isotropy, yields 
a scattering cross section which is not defined along the direction of 
primary propagation. Thus, their function cannot be used as the integrand 
in equation ( 6 ) . 
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E 



»y p4Jt 

/ au^oo dy 


(7) 


Because of assumptions evoked by Booker and Gordon in the derivation of 
equation ( 7 ) is rigorously valid only for cases where an element of 
thickness of the turbulent stratum dy exceeds the integral scale Z. 

In other words, the tTirbulence must be sensibly homogeneous throughout 
the volume (dy)^. Further, the total thickness of the turbulent layer, 
y, must exceed the radiation wave length, A. Finally, it should be noted 
that the effects of multiple scattering have been neglected. 


It is known (refs. 1^ and I 5 , for example) that photometric experi- 
ments involving Lambert^ s law can be subject to error because it is 
impossible to restrict the angular aperture of the measuring instrument 
to zero. Thus, an instrument centered on the axis of primary propagation 
and situated a finite distance from the scattering volume gathers, in 
addition to the attenuated plane wave, that portion of the scattered light 
lying within its field of view. To account for this effect, equation (7) 
is modified to read: 


where Ejyj is the intensity measured by an instrument having a field of 
view defined by the solid angle The instrument is considered located 

sufficiently far from the scattering layer that changes in solid angle 
with layer depth can be neglected. (This is the paraxial ray assumption 
of geometrical optics.) For an instrimient with an axially symmetric field 
of view alined with the axis of primary propagation, the solid angle sub- 
tended, (jj, is related to the scattering angle, 0, by 


Measurement of Light Intensity in Scattered Field 



y=5 4jt 


(8) 



(9) 
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After the foregoing equalities are used to effect substitutions in 
equation (5)> 'tlie angular integration indicated in the ejq)onent of equa- 
tion (8) can be carried out by elementary methods. The result is: 



where K = The scattering coefficient Ug (eq. (6)) corre- 

sponds to equation (lO) with = 0. 

In the case of velocity fluctuations (refs. l6, 17> l8) , it has 

been found that the integral scale is of the order of half the boundary- 
layer thickness. It is reasonable to expect that the scale of the density 
fluctuations will be of a similar order of magnitude, say about 1 inch, 
for flight and wind-tunnel boimdary layers . Because the wave length of 
yellow light is about 2x10“^ inch, the discussion can therefore be 
restricted to cases where I is very much greater than A. Consequently, 
only the leading term of equation (lO) contributes significantly to the 
scattering coefficient, and one obtains with the help of equation (9): 


a, 


. 2jt^Z /Ae^^ 


SM 


1 . 2 M 

1 - Sin^ 


1 + 


A 


sin 


2 ^ 


and, for 


Ug - 


2jt^Z 




( 11 ) 
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Furthermore, the results of references l6 and l8 indicate that the integral 
scale of the velocity fluctuations is, for practical purposes, independent 
of position in the boundary layer. When it is assumed that the integral 
scale of the density fluctuations is likewise independent of position, the 
integration with respect to distance indicated in equation (8) yields: 



and, for = 0: 







( 12 ) 


Zn 




J 


The first of equations (l2) relates the measured intensity ratio Ej^/Eq 
for given wave length of the incident radiation and structirre of the 
turbulence as characterized by integral scale and intensity of density 
fluctuations with the changes in angular aperture, 9^, of the viewing 
instrument. These results are applicable to turbulent compressible bound- 
ary layers subject to the following qualifications: (l) the integral scale 

is considered invariant throughout the layer; (2) the integral scale is 
very much greater than the radiation wave length; (3) the turbulence can 
be considered isotropic; and (4) the turbulence is characterized by a space 
correlation of the density fluctuations that obeys an exponential decay 
law. 


Photometric Exploration of Fluctuating Density Fields 


Although the foregoing analysis presimies that the- integral scale and 
intensity of the density fluctuations are known and the radiation intensity 
is the dependent variable, the converse is usually true in practice. It 
is therefore convenient also to regard equations (l2) as defining the 
turbulence structure in terms of known radiation variables , thereby per- 
mitting the extraction of useful information about tiorbulent layers. To 
this end the first of equations (12) can be rearranged to read: 
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It appears from equation ( 13 ) that photometric surveys of the light field 
after interaction with the tiirhulent layer, wherein the ratio E^/Eq is 
measured as the angular aperture fijyj of the photometer is varied, can 


yield values of the unknown turbulent parameters I and 


V=6 

I 


ciy. 


This is true because a graph of the left-hand side of equation ( 13 ) against 
the aperture function on the right-hand side yields determinate functions 
of the unknowns as the slope and intercept of a straight line (sketch (a)). 



Sketch (a) 

As part of the present investigation, the functional adequacy of equa- 
tions ( 12 ) and ( 13 ) when applied to turbulent compressible boundary layers 
was tested by comparison with photometric data. 


APPARATUS AND PROCEDURE 
Wind Tunnel 


The present experimental investigation was conducted in the Ames 1- 
by 3-foo't supersonic wind tunnel no. 1. The natural boundary layers that 
exist on the tunnel side walls were thickened artificially to obtain a 
range of test conditions. This closed circuit, continuous operation, 
variable pressirre wind tunnel is equipped with a flexible-plate nozzle 
that provides a range of supersonic Mach numbers from 1.20 to 2 . 5 O. A 
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range of subsonic Mach numbers from 0.40 to O.85 can also be provided at 
the test section by locating the minimiun area in the circuit downstream 
of the test section. 

The working fluid is dry, filtered, atmospheric air. Contamination 
within the circuit due to oil and manometer fluid vapors and leaks from 
the atmosphere is minimized by continuous air interchange at a rate suffi- 
cient to prevent water condensation in the nozzle. The absolute pressure 
in the settling chamber can be varied from I/5 atmosphere to 2 atmospheres 
to provide changes in Reynolds n\miber. Corresponding absolute tempera- 
tures vary from 520° R to 600° R and depend primarily upon the absolute 
pressure. Stagnation pressure and stagnation temperature can be measiired 
with a relative error of ±1 percent in the most uncertain cases. 

Some heat flow always takes place between the room and the nozzle 
boundary layers, but the heat-transfer rates are small. This is manifested 
by the fact that the wall temperatures have equilibriim values within a 
degree or two of the theoretical recovery temperatiires of the stream. The 
temperatures on the test -section wall vary from about 20° F above room 
temperature at the lowest Mach number and highest pressure to about 25° F 
below room temperature at the highest Mach number and lowest pressure. 
Although it is usually possible to select a stagnation pressure which pro- 
duces zero heat transfer at any given Mach number, it was felt that a close 
control of heat transfer was unwarranted because the heat flow rates are 
small . 


Instrumentation 


Aerodynamic .- Boimdary- layer pressure data were provided by either 
^ 28-tube or a 36-tube total -pressure rake which spanned the test section 
on the horizontal median plane, 4-5/8 inches downstream from the center 
of the 9-inch-diameter viewing windows (fig. 2), and also by 25 static- 
pressiire orifices on each test -section wall. Pressures were indicated by 
a conventional liquid manometer system. In addition, three hot-wire ane- 
mometer probes were installed at the rake station with the wires normal 
to the air stream. One probe was mounted integrally with the rake on the 
wind-tunnel axis. The remaining two were installed on a horizontal trav- 
ersing mechanism so that they could be traversed through the boundary 
layers on opposite walls of the wind tunnel. The hot-wire anemometer was 
used for obtaining relative boundary- layer and free-stream turbulence 
levels and for detecting the location of the edges of the boundary layers . 

The relative error in determining Mach number from the rake total 
head and wall static pressirre is ±1 percent. Locations of the centers 
of the 1/32 inch outside diameter total head tubes were measured to the 
nearest l/6k inch. 
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Optical .- The optical equipment used in the present investigation 
was substantially the same as that employed by Baskins and Hamilton 
(refs. 5 and 6). The principal components were three identical f/8.6 
Newtonian telescopes having focal lengths of 21.45 inches. Two of the 
telescopes were used as source collimators, and the remaining instrimient 
was equipped as a telephotometer. Light from battery-operated tungsten- 
filament incandescent bulbs was relayed to 0.00025- inch-diameter apertirres 
at the principal foci of each collimator . Electronic photometers viewing 
the sources directly served as monitors of the two source intensities. 

As is shown in figure 2 a collimator was arranged on each side of the 
wind tunnel to direct, as required, plane waves toward either of the 
high-quality windows (see ref. 6 for description) in the wind-tunnel side 
walls. A half-silvered mirror and the receiving telescope were situated 
on opposite sides of the wind tunnel so that some of the radiant energy 
from either collimator could enter the receiver after passing either once 
or twice through the air stream.^ After passing through the final wind- 
tunnel window, light from the selected collimator entered the receiver 
through a plane -parallel glass plate which capped the telescope barrel 
at an angle of 45° with respect to the optical axis. The reflex mirror 
which is characteristic of the Newtonian arrangement was situated at the 
center of this plate and consisted of an elliptical deposit of aluminum 
having a nominal minor diameter of l/4 inch. Thus the entrance pupil was 
an annular window having a nominal major diameter of 2-l/2 inches and a 
nominal minor diameter of l/4 inch, and the resulting diffraction pattern 
at the principal focus was of the center-blocked Airy type. Although the 
source diameter was not sufficiently small to meet the rather strict 
coherence criterion of Hopkins (ref. 19), it was about one-half the maxi- 
mum diameter considered adequate by Gabor (ref. 20) and two-thirds the 
diameter considered a good working compromise by Rogers (ref. 2l) . 

The receiving telescope was equipped with means for inserting at the 
principal focus a series of pinhole apertures having diameters ranging 
from 0.0007 to 0.0330 inch. When in position, the center of each pinhole 


^Because the receiver and source collimator on the same side of the 
tunnel were physically separated, the double-transmission beam incident 
on the mirror irradiated slightly different portions of the boundary layers 
before than it did after reflection toward the receiver. The maximum 
displacement between incident and reflected beams was about 1 inch at the 
STorface of the near window. 

The original optical arrangement for the double-transmission tests 
employed an opaque first -surface mirror which could be inserted at the 
location of the half- silvered mirror. Because installation and calibra- 
tion became a lengthy and frequently repeated process, and because the 
wind tunnel had to be shut down during changes from single-transmission 
to double -transmission operation, the arrangement using the half- silvered 
mirror proved to be the more desirable from an operational standpoint . 
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coincided with the optical axis of the receiver. The radiation which 
passed each pinhole was relayed by a lens system as a collimated beam 
to a l/8-inch-diameter spot on the cathode of a 1P21 photomultiplier tube. 
The photomultiplier tube acted as the transducer of a commercially avail- 
able indicating photometer. Repeated calibrations on an optical bench 
showed that this instrument was stable with time and had a linear response 
within about O .5 percent of full scale for output readings below I 5 micro- 
amperes. The combination of light source, optical components, and photo- 
electric cell possessed a spectral sensitivity having a maximum response 
at a wave length of about 5200 angstrom iinits. The energy intercepted by 
the phototube was measured with a maximimi relative error of ±2 percent. 

With the various pinholes in place, the angular aperture of the 
photocell varied from 1.7X10“^ to radian. Aperture diameters 

were measured to the nearest 0.00001 inch, which corresponds to slightly 
less than 0.03X10”^ radian in angular aperture. Angular apertures cor- 
responding to the various pinholes are recorded in table I. 

The collimators and telephotometer were mounted on a cradle which 
was isolated from the floor and the tunnel structure by an elastic cord 
suspension. The cradle assembly and test section were housed in a cabin 
to minimize extraneous illumination and room turbulence. 


Tests 


Aerodynamic parameters .- In addition to the individual effects of 
Mach number and free- stream density on the diffusion of light, the effect 
of bound ary- layer thickness was determined. The 1-1/2- to 2-inch-thick 
natural boundary layers on the wind-tunnel side walls were artificially 
thickened by two methods. The first was a boundary- layer trip consisting 
of a series of adjustable-height fences transverse to the side-wall 
boundary- layer flow in the low subsonic region at the nozzle entrance. 

At subsonic Mach numbers this arrangement caused thickening by factors 
up to 3* However, the effectiveness of the trip as a thickening device 
fell off rapidly with increasing Mach number, and practically no thicken- 
ing occurred at a Mach number of 2.5 • The second method consisted of 
applying a distributed roughness coating that extended in ^-±nch wide 
bands on each wall to various lengths down the nozzle. The roughness 
consisted of irregular- shaped grains (coarsely ground coffee beans having 
maximum dimensions ranging from I /16 to about I /8 inch) attached to the 
steel walls with a rubber -base adhesive. With the roughness bands extend- 
ing to the viewing windows, the boundary layers at a Mach niimber of 2.5 
were about 3 -l /2 inches thick, as opposed to the natural thickness of 
1 - 3/4 inches. However, at subsonic Mach numbers the distributed rough- 
ness was not so effective in increasing total thickness as the boundary- 
layer trips. 
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Some observations were made through separated boundary layers. With 
the M = 1.75 nozzle setting, driving the wind-tunnel terminal shock wave 
structure into the test section caused the boundary layer to separate. 

The shock was positioned by adjusting the back pressure in the diffuser 
so that the region of maximum separation (l/2 to 1 inch thick) was over 
the observing windows. 

Although relative comparisons with the hot-wire anemometer indicated 
that the free- stream fluctuation level at all Mach mmibers was at least an 
order of magnitude less than that in the boundary layers, there existed 
the possibility that the free-stream fluctuations could contribute measur- 
ably to the optical diffusion because, with thin boundary layers, the 
optical path through the free stream was up to four times as long as the 
path through the boundary layers themselves. This possibility was explored 
by employing a turbulence -promoting grid of bars at the nozzle entrance 
to increase the free-stream fluctuation level. Approximately doubling the 
free-stream fluctuation level as measured by the hot-wire anemometer had 
a negligible effect on the optical data at Mach numbers of 0.8 and 2.5* 

A limited number of tests were conducted to determine the effect of 
air-stream humidity on the optical transmission characteristics. At super- 
sonic Mach numbers, so long as the water vapor did not condense, little 
change in image quality could be detected. However, when condensate foimied 
in the test section, the effect on the image was as if the source intensity 
had been decreased. This reduction of intensity due to water vapor, which 
amounted to as much as 20 percent, suggests that processes of skin cooling 
by liquid evaporation must be carefully controlled if liquids are injected 
upstream of receivers of visible radiation. Nominal test conditions at 
which the bulk of the data was obtained are given in table II. With few 
exceptions, any actual test condition corresponded to one of these nominal 
conditions with a relative error of less than ±3 percent. 

Optical parameters *- To avoid the experimental difficulties associated 
with establishing a procedure for viewing through a single turbulent bound- 
ary layer, the two boundary layers on the wind-tunnel side walls were used 
as the basic configuration. It was felt that this experimental simplifi- 
cation was justified because equation (2) gives the possibility for reduc- 
ing the results for two boundary layers to an equivalent result for one. 
Moreover, the validity of equation (2) is amenable to experimental verifi- 
cation by comparison of the results of observations through two boundary 
layers with the results of observations through four (eq. (3))* Thus, 
for each setting of Mach number, density level, and boundary- layer thick- 
ness, sets of optical data were taken for the two cases of single and 
double transmission of light beams through the test section. 

A set of optical data consisted of the series of measurements of 
radiant flux which crossed the focal plane of the receiving telescope 
through each of the pinhole apertures listed in table I. The output 
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indication of the receiving photocell was directly proportional to the 
time-averaged value of the radiant energy which passed through each pin- 
hole per imit time, and was related to the time-averaged value of the 
local intensity as follows: 


P 


F 




(lU) 


where d is the pinhole diameter, Ip is the time-averaged local inten- 
sity, considered a function of radius alone, and Ip is the average 
intensity over the aperture area. 


To minimize the effects of molecular scattering, and imperfections 
in the optical components, and, further, to account for diffraction in 
the telephotometer, equation (lU) was normalized with respect to condi- 
tions for no flow in the wind tunnel (test section at atmospheric pressure 
and temperature) : 


Pp Ep Ip 

%F %F ^NF 
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(15) 


According to equation (15) when Ip is a different fxmction of r than 
I]^p, the ratio Ep/Ejjp becomes a function of the pinhole radius, d/ 2 , 
and thus of the angular aperture since 
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2 F 


_d_ 

2 F 


(16) 


where F is the focal length. As is shown in reference 6 , equation (15) 
effectively eliminates the diffraction effects brought about by the finite 
entrance pupil of the telephotometer, and thus the quantity Ep/Ejjp is 
a measure of the flux within a given solid angle due to boimdary- layer 
effects alone, referred to the total flux per unit area, Eq^ in the 
entering wave front; that is. 
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Data to check equation (I7) were obtained by varying the focal ratio of 
the telephotometer from f/8.6 to f/86. Decreasing the entrance pupil 
diameter from 2.5 inches to 0.25 inch increased the size of the diffraction 
image by factors up to 10. Tests with the largest focal ratio were incon- 
clusive because the relatively small quantity of light available was 
spread over a large area of the focal plane, and difficulties were enco\m- 
tered in accurately centering the pinhole apertures. However, changes in 
focal ratio up to f/'^Z had no measurable effect on the ratio 
for given flow conditions. 

Brief tests to assess the dispersion according to wave length in the 
scattered fields (Tyndall spectra) by capping the receiving telescope with 
filters were inconclusive because the spectral range covered by the avail- 
able filters was small. It should be pointed out, however, that Baskins 
and Hamilton (ref. 6) report definite evidence of a spread in the scattered 
fields according to wave length, the shorter wave lengths being dispersed 
through larger angles. 

Although the bulk of the data was obtained at normal incidence (cp = O) , 
limited tests were also carried out at an incidence angle of i^5°. Because 
of inadequate clearances between the tunnel structure and the rearranged 
optical components, data could be obtained only for the case of double 
transmission. 

Although the image of the point source was, typically, symmetric 
about the optical axis, a number of extraneous aerodynamic effects were 
encountered which had an undesirable influence on the image symmetry . 

If the incident beam intercepted a weak shock wave, the image acquired a 
tail oriented at the shock angle. Such situations were remedied by modi- 
fication of the nozzle contours to eliminate the disturbance. During 
initial testing near a Mach number of O.65, the image was observed to 
break up into a coarse -grating- type diffraction pattern. This behavior 
was traced to an acoustic grid of ultrasonic standing waves in the wind- 
tunnel test section excited by the boundary- layer rake. When the rake, 
which had a bl\mt trailing edge, was modified by addition of a sharp 
trailing edge, the image again became axisymmetric . 

The fact that the test- section walls could not be maintained at room 
temperature for all test conditions gave rise to free-convective flows 
over the exteriors of the viewing windows. It is believed that in most 
cases this flow was laminar, but there was no convenient way to effect a 
calibration and it is conceivable that in certain instances free turbu- 
lence over the windows might have made a small contribution to the image 
deterioration. Accumulations of dirt on the windows were a continuing 
source of difficulty and annoyance. However, calibrations through clean 
and dirty windows established that runs of U to 6 hours duration intro- 
duced additional losses of less than 1 percent. 
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The optical parameters which were varied in combination with the 
various aerodynamic conditions are given in table III , and table IV is 
a tabulation of the aerodynamic and optical data for all configurations. 


RESULTS AIv'D DISCUSSION 
Photometric Measurements 


Effect of flow on transmitted power .- With the largest pinhole aper- 
ture (d = 0.033 in.) in the focal plane of the receiving telescope it was 
found that there was little difference in the radiant power received by 
the photocell for a condition of no flow in the wind t\mnel and for the 
flow condition corresponding to the greatest image deterioration. In 
other words, the turbulent contributions to absorption, side- and back- 
scattering were not measurable with the present instrumentation, and the 
total energy in the incident wave for practical purposes always passed 
the entrance pupil of the telephotometer irrespective of the flow condi- 
tion in the wind tunnel. This result, which was first reported by Baskins 
and Hamilton (ref. 6 ) can be interpreted to indicate (eq.(5)) that, in 
boundary layers having total thicknesses of the order of 1 inch, the 
integral scale of the density fluctuations is very much greater than the 
effective radiation wave length of 5200 angstrom units. Thus, in accord- 
ance with the theory of Booker and Gordon, it can be stated that the light 
field scattered by a compressible turbulent boundary layer is beamed prin- 
cipally in the forward direction. The maximum scattering angle encoun- 
tered in the present tests has a value of O.OOO 6 radian, which corresponds 
to a condition of no loss through a pinhole aperture having a diameter of 
0.025 inch. 

Comparison of one- and two-transmission data .- In figure 3 the per- 
centages of light reaching the phototube through various pinholes'^ for 
various flow conditions after penetrating two boimdary layers are plotted 
against the corresponding percentages received after penetrating four. 

Only the data taken with the half- silvered mirror in place are shown. 

All these data, irrespective of flow condition or angular aperture, cluster 
with remarkably small scatter about a single ciirve. However, the mean line 
passing through the experimental data differs systematically from the theo- 
retical prediction based on Lambert's law (eq. (3))* From equations (12) 
it is apparent that Lambert’s law should apply for given angular aperture, 
because doubling the number of identical boundary layers penetrated by 
the light beam merely doubles the value of the fluctuation integral, 
thereby squaring the original intensity ratio as was sho^vn previously 
(eq. ( 3 )). Although part of the discrepancy between experiment and theory 
might be attrib utable to the fact that the paths of the incident and 

■^Because of negligible losses, the data for the two largest pinholes 
have been omitted from figure 3 "the following figures. 
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reflected light for the double-transmission cases did not coincide, the 
discrepancy is believed more probably to be caused by two secondary 
effects that are neglected in the theoretical treatment, namely, second- 
ary scattering and scattering from rough interfaces between the boundary 
layers and the free stream. In the region where from 10 to 65 percent 
of the total energy is received through two boundary layers, a maximum 
increment of 0.04 in the value of (Ep/Ej^Tp)^ over corresponding predictions 
based on Lambert’s law occurs. This increment is believed due to recep- 
tion of light which has been scattered more than once, that is, from solid 
angles outside that of the photometer. Further, in the region where from 
65 to 100 percent is received, a maximum decrement between experiment and 
theory of 0.02 occurs. This decrement is attributed to a scattering of 
energy from the rough interfaces which occur between the boundary layers 
and the free stream (ref. 17)* Due to interface scattering alone, the 
experimental data would be expected to lie below the theoretical curve 
throughout the range of intensity ratios and to be a function of the flow 
conditions. Secondary scattering, however, appears to act in a direction 
to counteract the decrement due to interface scattering and, moreover, 
to exert an influence which increases as the intensity ratio decreases. 

The net result of these two secondary effects is to cause the data of 
figure 3 to be independent of angular aperture and flow condition, but 
to depart systematically from the theoretical curve. 

Because of these secondary effects, an apparent change in integral 
scale, Z, will be observed, and the corresponding logarithms of the 
intensities at = 0 v/ill not, in general, be in the ratio of 2 to 1 
for data corresponding to the single- and double-transmission cases. 

Thus after equation ( 13 ) is written for two and four boundary layers and 
is eliminated between them, an expression is obtained that contains 
the apparent integral scales and limiting intensities from the two 
measurements : 


1 




(18) 


where 
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The empirical c\irve defined by the data in figure 3 is given by: 



E 


0»44— 0»14 Zn(Ej^/Eo)2 


0/2 


for 



> 0.10 


Therefore , the data indicate the following empirical result 


1 + K A Zn(E/Eo)2 
1 Zn(E/Eo)4 


O.ll = i A + K4 _ 

Zn(E/Eo)4 + ^2 

Consequently, the integral scale and intensity of turbulent density fluc- 
tuations for two boundary layers can be obtained in terms of the corre- 
sponding quantities for four boundary layers for this series of experiments 
by use of the following empirical expressions (for (E}^/Eo )2 > 0.10): 



Zn 



0.U4 Zn(E/Eo)^ 

1 + 0.14 Zn(E/Eo)4 


1 + 



1 + (4itZ4A)^ 

1 + 0.14 Zn(E/Eo)4 






( 19 ) 


It is reasonable to expect that for the same intensity ratios, a 
comparison between one boimdary layer and two boundary layers would yield 
smaller departures from Lambert »s law than is shown in figure 3 for the 
case of two and four because the mmiber of boundary- layer interfaces is 
cut in half and because opportunity for the occurrence of secondary scat- 
tering decreases. The dash-dot curve lying halfway between the experi- 
mental curve and the prediction based on Lambert’s law, which is shown in 
figure 3^ is an estimate of the relation expected to hold between one 
boundary layer and two. This curve is subsequently used to reduce the 
optical data to results applicable to an equivalent single boundary layer. 
The empirical equation for this curve is: 
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0.48-0.06 Zn(EM/Eo)i 
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(%) > 0.28 


( 20 ) 


Comparison of data with scattering theory .- Figure has been pre- 
pared to compare the results of the theoretical scattering analysis with 
the test results for compressible turbulent boimdary layers. The theory 
based on the two-parameter turbulence model of Booker and Gordon predicts 
the straight-line functional dependence given by equation (l 3 ) and depicted 
in sketch (a). Figure 4(a) presents single -transmission data at a Mach 
number of 2.00 and a density ratio of 0.46 for each of boundary- layer 
configurations 0 through (0 (table II ) plotted in the coordinate system 
suggested by the analysis. In each case the straight-line dependence 
predicted by theory is obtained, and thus the turbulent density fluctua- 
tions can be defined in terms of the two parameters given by the extra- 
polateci intercepts with the axes. 

That wide changes in Mach niJinber and density ratio have little effect 
on this fmctional behavior is further illustrated in figure 4(b) , which 
is a normalized form of figure 4(a). The ordinate of figure 4(b) is the 
left-hand side, and the abscissa is the right-hand side of the following 
rearranged form of equation (l 3 )* 


T 1 1 + (4jtZ2A)^ 2 % 

1 1 = tan —— 

Zn(%/Eo)2 Zn(E/Eo)2 ^n(E/Eo)2 ^ 


Optical data corresponding to a given flow condition were plotted 
cated in figure ^(u) and fitted with a straight line. The x and 

intercepts determine the constants Zn(E/Eo )2 jn(E/E ) 


as indi- 
y axis 

Graphi- 


cal curve fitting was used for most cases; hut where the scatter was 
considered excessive, the constants were determined by the method of 
least squares. Fitting the straight line was accomplished by the method 
of weighting of functions described in reference 22. In the present case 
this weighting procedure had the effect of decreasing the weight of those 
data corresponding to intensity ratios which were both larger and smaller 
than a value of l/e (= 0 . 368 ). 


Notwithstanding the considerable variation of Mach number, free- 
stream density, and boundary- layer thickness represented by the various 
symbols shown in figure 4, it can be observed that departures from line- 
arity are small. It therefore appears that the fluctuating density fields 
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in compressible tiirbulent boundary layers behave as though they were iso- 
tropic insofajT as the transmission of light is concerned. Consequently, 
photometric measurements in the radiation field produced by the inter- 
action of a plane light wave with a tirrbulent boundary layer can be used 
in conjunction with scattering theory to deduce average values of the 
integral scale and the intensity of turbulent density fluctuations. 


Effect of Turbulent Density Fluctuations on Optical Images 


It has been shown so far that the scattering analysis based on the 
cross section of Booker and Gordon is in good fimctional agreement with 
experiment, and, further, that the turbulence structure can be defined 
by two parameters deduced from the data. With this information, one can 
compute for a given flow the distribution of radiant intensity at the 
focal plane of a diffracting optical imaging device. To carry out this 
calculation, one equates equations (l2) and ( 15 ) by using equation (l?)^ 
and changes variables as follows: 





After differentiating the resulting expression with respect to the dimen- 
sionless radius | and simplifying, one obtains 


Zn e/Eq 


= 
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For an axially symmetric objective one has for Inf well-known result 

of Airy: 


Inf 


Inf 


5=0 


'2Ji(5) ~ 
_ I _ 


where Ji(|) is the Bessel function of unit order. It therefore follows 
that the intensity at a point | in the focal plane due to turbulence 
relative to that on the optical axis without turbulence is: 
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Graphs of equation (2l) for the no-flow condition and for three 
typical flow conditions with boundary-layer configuration are compared 
in figure 5 with each other and also with corresponding photomicrographs 
of the single-transmission image at the focal plane of the receiving 
telescope. For each graphs the two parameters involving the turbulence 
structure were evaluated from single-transmission photometric measurements 
in the manner previously described^ and the two optical constants, A and 
D, were assigned the values appropriate to the instrumentation used in 
the present tests. The effect of center-blocking is small and can be 
safely neglected. 

The ordinates at the point ^ = 0 in each of the four graphs in 
figure 5 have been given the same physical height. However, the relative 
n-umerical values of the peaks are actually shown to be in the ratios 
1:0.5:0. 2 :0.1 as one proceeds down the col-umn. An intensity ratio of 
unity on the optical axis corresponds to the no-flow condition shown by 
the graph and the photographs along the top row. Each of the graphs of 
the intensity ratios for the three flow conditions also shows as an addi- 
tional dashed curve an attenuated no-flow distribution. This dashed 
curve arises from the first term on the right-hand side of equation (2l) 
and can be interpreted as the portion of the intensity distribution com- 
posed of unscattered radiation. Note on the graphs that as the value of 
the central maximum decreases with increasing air density, the height of 
the total intensity distribution curve relative to the corresponding no- 
flow curve increases and that the dark rings of the diffraction pattern 
(l = 3-8^ 7*0, 10.1, . . .) disappear. This indicates that increasing 
fractions of radiant energy are scattered outside the Airy disc. For 
example, the ratio of the flow and no-flow ordinates at the maximum of 
the second bright ring (^ = 8.6) are l/l, 22/l, and 32/l^ respectively, 
for the 0.157^ 0*295^ 0-379 density ratios. 

The photomicrographs of the images to the right of each plot in 
figure 5 sire visual evidence of the increased scatter that occurred as 
the density was increased at a Mach number of 2. Unfortunately, due to 
halation in the film and limitations in the reproduction process, the 
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photographic images show considerably less detail than can be observed 
visually^ or for that matter^ on the original negatives. The Airy disc, 
the first two bright rings, and portions of the third bright ring can be 
readily distinguished on the negatives for the case of no flow. The 
extraneous spots of light visible on some of the photographs are "ghosts” 
originating from second- surface reflections, and they represent less than 
1 percent of the total flux. The exposure times (l5^ ^5^ and 90 seconds) 
for the three columns of images shown were selected so that points of 
equal brightness on the film in any row were exposed by light intensities 
in the ratio 6:2:1. With flow at a Mach number of 2 (bottom three rows 
of pictures), the intensity of the central disc is noticeably attenuated 
(15-second column) and the zone of the focal plane outside the disc becomes 
illuminated correspondingly (45- and 90-second coliomns) . Except for the 
small influence of secondary scattering, these qualitative effects of the 
scatter and redistribution have been predicted by the previous analysis. 

A further feature of the photographs for the flow conditions shown 
in figure 5 should be mentioned. The maximum intensity ratios calcu- 
lated as belonging to each row and the selected exposure times for each 
column of photographs cause adjacent pictures on a diagonal (dhl, for 
example) to exhibit approximately the same response; for by the reciproc- 
ity law for photographic emulsions, the increased exposure times along 
the diagonal serve to counteract the decrease in light intensity caused 
by the additional scattering that occurred. Therefore, increases of the 
diameter of the image along the diagonal dhl illustrate that even 
though one can compensate for the attenuation due to turbulence by 
increasing the gain of the system (that is, increasing exposure time here 
or increasing the amplification factor of a photomultiplier circuit) a 
loss in resolution will be inevitable. 


Dependence of Optical Transmission Characteristics on 
Mean Flow Conditions 


Establishment of correlating parameter .- There has been no attempt 
so far to connect the scattering of radiant energy with the dynamics of 
the mean boundary- layer flow because it is clear from the preceding dis- 
cussion that Mach number, density level, and boundary- layer thickness are 
important only insofar as they contribute to the formation of density 
"lumps.” For reasons of practical utility, however, it would be desirable 
to utilize some average characteristic of the mean boundary- layer flow 
as a measure of the turbulence pattern controlling the scattering of light. 
Unfortimately, the lack of knowledge concerning the relationship between 
the turbulence structure and the mean flow in a compressible turbulent 
boimdary layer prohibits the formulation on any other than an empirical 
basis . 
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As a preliminary to the formulation of such a parameter, it is noted 
that, when plotted as shown in sketch (a) and figure 4(a), the photometric 
data exhibit the following characteristics: 

1. Increases of Mach number, density, and bo\mdary- layer thickness 
produce substantial decreases in the slopes of the straight lines passed 
through the data. 

2. Relatively large changes in Mach number, density, and boundary- 
layer thickness produce comparatively small changes in the absolute values 
of the extrapolated intercepts with the abscissa, and, moreover, the 
intercepts in most cases have small numerical values as compared to the 
nijmerical values of the apertxire function tan (6j^/2). 

In other words, displacements of the present data regarded as shifts 
in the coordinate axes are much less extensive than displacements regarded 
as rotations of the coordinate axes. Thus the slopes of the lines alone 
will serve as a good over-all measirre of the effects of Mach number, 
density, and boimdary- layer thickness on the scattering process, and one 
can, with little loss in accuracy, consider that the straight lines passed 
through the data also pass through the origin of the coordinate system. 

The mathematical conditions to be met for this simplification to be valid 
are: 


In the present tests these limitations are fulfilled for all of the pin- 
hole apertures for some of the test conditions, and for most of the larger 
pinhole apertures for all test conditions. Equation (l2) can therefore 
be written in the following form independent of explicit dependence on 
the radiation wave length: 



or 




6 


( 22 ) 


Thus, a parameter based on mean flow conditions which is proportional to 
the right-hand side of equation (22) can be expected to effect a correla- 
tion of the intensity ratio data, at least for the larger pinhole 
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apertures. It should he noted that the flow-dependent portion of the 
exponent of equation (22) can he cast into a form which hears a very- 
close resemblance to the corresponding expression obtained by Liepmann 
(ref. 2) for the mean square deviation of a ray. If one rewrites 
equation (22) as follows 



and notes that the left-hand side must be constant for a given flow con- 
dition, it is obvious that each set of optical data can be described by 
a single number, say the mean square deviation of a ray, 0^. If, fiirther, 
the integral on the right-hand side of the above expression is written 
as the product of an average value of the integrand and the boimdary- layer 
thickness, one obtains 


0 ^ ~ 


In 




(23) 


This relation corresponds f-unctionally to equation III-2 of reference 2. 
Liepmann approximated equation (23) in terms of the mean flow variables 
by: (1) setting the root mean square density fluctuation proportional 

to the over-all change in mean density across the bo\mdary layer, 

(2) setting the scale of the fluctuations proportional to the thickness 
of the laminar sublayer, and (3) setting the total boundary- layer thick- 
ness proportional to the length of boundary- layer run divided by the I/5 
power of the length Reynolds number. This approach leads to the following 
result for the root mean square deviation of a ray: 




' (7-1)Mqo72 ~ 

.1+(7-1)Moo72- 


/ length 
V^Reynolds no 


local skin-friction'\^^ ^ 
coefficient J 


For any one of the boundary- layer configurations listed in table II the 
optical data from the present experiments can be correlated by a parameter 
proportional only to the first two factors in the foregoing expression, 
that is, the free-stream density and free-stream Mach number function. 
However, attempts to incorporate into the parameter the relatively weak 
effects of the various boimdary- layer thicknesses and profile shapes 
(Reynolds number and skin-friction coefficient) by the use of functions 
of displacement, momentum, or energy-loss thickness were not an unquali- 
fied success. These correlations usually led to a spread of the 
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imseparated-flow data as a fimction of Mach number and did not correlate 
the separated boundary- layer data at all. 

The inadequacy of that portion of Liepmann’ s parameter proportional 
to the ratio of bomdary-layer thickness and turbulence scale is believed 
due to the use in equation (23) of flow parameters proportional to the 
microscale of the tiorbulence. The present analysis, however, identifies 
the length, I, in equation (23) as the integral scale of the fluctuations 
rather than as the microscale. According to references I6, I7, and I8, 
the flow parameter proportional to the integral scale is the total 
boundary- layer thickness, &. Hence the right-hand side of equation (23) 
is actually insensitive to bo\mdary- layer thickness except as it influ- 
ences the local intensity of the density fluctuations and. thus the aver- 
age. This latter effect can be taken into accoimt by a slight generali- 
zation of that portion of Liepmann' s parameter which includes the fluc- 
tuation contributions. Instead of setting the fluctuations proportional 
to the over-all difference between the density at the wall and in the 
free stream as Liepmann has done, it has been fo-und advantageous to employ 
the Integral across the layer of the difference between the free-stream 
density and the local density within the layer. Thus, defining the 
boundary- layer parameter: 



and introducing the usual flat-plate assumptions that the total tempera- 
tirre and static pressure do not vary across the boundary layer, one finds 
that 
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This relation is rendered dimensionless by introducing the standard sea- 
level density, p, and the effective radiation wave length. A, so that 
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With respect to the free-stream density and free-stream Mach nimber 
dependences, the parameter p (eq. ( 24 )) is identical to the analogous 
expression given by Liepmann in reference. 2. However, whereas Liepmann' s 
parameter incorporates the ratio of the total boundary- layer thickness 
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to turbulence scale in terms of Reynolds number and skin-friction coef- 
ficient, the parameter (3 utilizes a thickness weighted with respect to 
the Mach number profile. No general significance can be attached to the 
parameter p, and its use can be justified only on the practical groimd 
that it is more effective in correlating the data than the other arrange- 
ments which were tried. 

Samples of the boundary- layer Mach niomber profiles typical of those 
used to calculate the parameter p are shown in figure 6, for a free- 
stream density of 0.23 standard sea-level atmospheres. The small dif- 
ferences between boundary- layer Mach numbers at corresponding y dis- 
tances on opposite walls of the wind tunnel were averaged arithmetically. 
Distances through the boundary layer are normalized with respect to the 
energy-loss thickness, which is defined in the figure. A sample of the 
Mach number profile for the case of the separated boundary layer (con- 
figuration ( 5 ) ) is also shown. It can be seen from the figure that the 
shape of the separated layer is significantly different from the others. 

It will appear shortly that this difference had an effect on the optical 
transmission characteristics which was adequately accounted for by the 
parameter p . 

Correlation of optical data with aerodynamic conditions .- Figure 7 
shows the fraction of energy remaining in the light beam after traversing 
four and two boundary layers as a function of the flow and optical vari- 
ables combined into the parameter p/Oj^cos cp. The limits for one boundary 
layer were estimated by the use of figure 3- Tbe correlation parameter 
incorporates the factor cos cp so that data obtained at the angle of 
incidence of 45^ could be included in the figure. 

Notwithstanding the wide range of Mach numbers, density levels, 
boundary- layer shapes and thicknesses represented by the symbols defining 
the family of curves in figure 7^ irrespective of the small departures 
from insulated wall conditions due to heat transfer, the parameter 
p/%cos cp effects an adequate correlation of all the data. The fact 
that this parameter correlates the data as shown implies that 



This proportionality follows from the two observations that the data in 
figure 7 can be empirically represented by an equation of the form 
Zn Ej^/Eq - - (p/0j^cos cp)^ and that equation (22) gives the corresponding 

analytic prediction. Thus these data serve to provide an approximate 
over-all relationship between the turbulence structure and the mean flow 
conditions in the boundary layers. 


30 


NACA RM A56B21 


Some of the data shown in figure 7 which are associated with the 
smaller pinhole apertures show a tendency to break off from the main group 
at various positions along the descending portion of the hand. Curves 
for a typical optical set (filled symbols) are shown illustrating this 
tendency. The dashed portions of the curves correspond to a linear extra- 
polation of the data; in figure 4(a) this region of extrapolation corre- 
sponds to the distance between the y intercept and the first data point. 
This behavior is in accord with the analytic predictions, which indicate 
that for constant turbulence structure, that is, constant p, a decrease 
of the angular aperture to zero causes the intensity ratio to approach a 
finite lower limit (fig. 4(a)). The exact manner in which the curve would 
approach the asymptote, however, is a fimction of the integral scale of 
the turbulence, and thus of the given flow condition. 

It is apparent from examination of figure 7 that the moderately thick 
turbulent bovindary layers explored in the present investigation, even at 
high subsonic Mach numbers, have an important influence on the propagation 
of light, the influence being to scatter increasing fractions of the origi- 
nally collimated beam as the Mach number, density, and boundary- layer 
thickness increase. In fact, the estimated c\irves for one boundary layer 
indicate that when p/0j^cos cp is about 5^10®, an attenuation of the inci- 
dent wave to a value of %/Eq = l/e ( =0.368) can be expected. In con- 
trast, the same attenuation in the sea— level atmosphere due to molecular 
scattering alone would require propagation over a distance of about 75 
miles. 

Effect of scattered light on optical resolution .- Although the 
tiorbulence- scattered light is deviated through very small angles, the 
deleterious effect on an optical imaging device can be large because 
scattered light from a given point in the object goes into a finite area 
in the image. The effect can be conveniently evaluated by establishing 
a resolution criterion analogous to the Rayleigh limit of resolution due 
to diffraction. 

The Rayleigh resolution limit (ref. 7) can be expressed as the angle 
subtended between two equally bright point objects whose images overlap 
in the back focal plane of an optical imaging instrimient by an amount just 
equal to the radius of the Airy disc of either image . Roughly 85 percent 
of the energy radiated into the entrance pupil by each point object is 
fo\md in each disc; and the intensity distributions and disc radii are 
fixed by the focal ratio of the instrument and the wave length of the 
radiation. The angle subtended at the objective of an axially symmetric 
instrument by the Airy disc radius is (for D >>A): 

0R . 


D 


( 26 ) 
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On the other hand, light which has passed through a turbulent boundary- 
layer has been seen to appear at the focal plane in a much larger disc 
(fig. 5 ) } with energy distribution determined primarily by the intensity 
and scale of the boundary- layer density fluctuations. In analogy to the 
Rayleigh criterion, one can propose that the radius of the disc containing 
85 percent of the total energy establishes the resolution limit due to 
turbulence. Thus, one accepts as tolerable a 15 -percent loss, which 
serves to identify in figure 7 a minim\mi value of p/ 0 jj^cos cp of about 
10® for a single boimdary layer. The ratio of the boundary- layer reso- 
lution limit for any value of p and the Rayleigh resolution limit due 
to diffraction for an axially symmetric imaging device then can be 
expressed as 


®R 


PP 

1.22A cos cp 


XlO 


(27) 


The interpretation of equation ( 27 ) is best illustrated by an example. 

The resolution limit of the 2.5-inch-diameter objective used in the present 
experiments is about 2 seconds of arc; that is, 1.22 A/d is about 1 x 10 “^ 
radian for A = 5200 angstroms. With p = 3900, the ratio of the resolu- 
tion limit due to turbulence to the resolution limit due to diffraction 
takes a value of 3.9. This value of p happens to correspond to a 
1-3/4- inch-thick turbulent boimdary layer at a Mach mmiber of 2.5 and a 
density altitude of about 4-5,000 feet (table IV). It should be noted that 
according to equation ( 27 ) an objective of large diameter will be more 
severely affected than will one of moderate size. Optical imaging equip- 
ment such as sextants, cameras, and fire control sights will be similarly 
affected. 


CONCLUSIONS 


The results of tests to determine the optical transmission character- 
istics of turbulent boimdary layers on an effectively insulated flat plate 
over ranges of Mach number from 0.4 to 2.5, free-stream density from 0.12 
to 0.93 standard sea-level atmospheres, and boundary- layer thickness from 
I-I /2 to 4-1/2 inches lead to the following conclusions: 

1. The loss in radiant intensity from a collimated beam of impolar- 
ized white light after penetrating a turbulent boundary layer depends 
mainly upon the integral scale and intensity of the density fluctuations. 
This result is in agreement with a theoretical prediction based on the 
scattering cross section of Booker and Gordon. A parameter proportional 


32 


NACA EM A56B21 


to the space integral of the cumulative change in mean density across 
the boundary layer is fo\md to provide an over-all measure of the inten- 
sity loss over the range of flow parameters investigated. 

2. The scattered light is deviated in all cases through small angles 
the maximum deviation from the direction of primary propagation being 
about 0.0006 radian in the present tests; and the angular distribution of 
scattered energy depends upon the integral scale of the density fluctua- 
tions. This result also agrees with the scattering theory based on the 
scattering cross section of Booker and Gordon. 

3. Comparison of the transmission characteristics of two boimdary 
layers with those of fo\rr boxondary layers indicates small deviations from 
Lambert's law of exponential attenuation which probably result from second 
ary scattering and diffuse refraction from the boundary- layer free-stream 
Interfaces. 

i<-. Due to the foregoing effects, significant losses in resolving 
power can be sustained by optical imaging devices which receive energy 
through compressible turbulent botindary layers. For example, a 1-3/1-- 
inch-thick turbulent bomdary layer at a Mach nimiber of 2.5 and a density 
altitude of 1-5,000 feet introduces a loss of resolving power greater by 
a factor of 3.9 than the Rayleigh resolution limit for a 2-l/2-inch 
objective. 

5. Photometric measurements in the radiation field produced by the 
interaction of a turbulent boundary layer with a plane light wave can be 
used in conjunction with electromagnetic theory to deduce average A/alues 
of the integral scale and the intensity of the turbulent density fluctua- 
tions. Provided that the limitations of the Booker and Gordon analysis 
are not exceeded, this conclusion is believed to be valid also for turbu- 
lent flows other than boimdary layers and for radiations at wave lengths 
outside the visible range. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Feb. 21, 1956 
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TABLE I.- RANGE OF ANGULAR APERTURES OF TELEPHOTOMETER 


Pinhole 

Diameter, 
dxlO®, 
in . 

Nominal angular 
aperture, 

radian 

Plotting 

symbol 

la, lb 

0.75, 0.71 

1.70 

0 

2a, 2b 

1.01, 1.07 

2.43 

□ 

3a, 3L 

1.28, 1.23 

2.91 

0 

ha, 4b 

1.80, 1.85 

4.24 

i> 

5a, 5L 

3.32, 3.27 

7.70 


6a, 6b 

4.20, 4.07 

9.63 

u 

7a, 7b 

4.86, 4.71 

11.13 

D 

8a, 8b 

8.09, 8.14 

18.92 

0 

9 

15.4 

35.9 


10 

33.0 

76.9 



TABLE II.- RANGE OF AERODYNAMIC VARIABLES 


Mach nc 
Ma 

). designation 
), Nominal 

I 

o.4o 

II 

0.67 

III 

0.80 

IV 

1.30 

V 

1.40 

VI 

1.75 

VII 

2.00 

VIII 

2.50 

Density 


Plotting 

0 

D 

C7 

17 

D> 

0 

□ 

0 

ratio 

p 

symbol 

(a) 

(a) 

(a) 

(a) 

(a) 

(a) 

(a) 

(a) 

A 

0.116 

-0 

— 

— 

— 

— 

©,© 

Q-0 

®-© 

Q<3 

B 

.232 

9 

©,(D 

©<D 

©-© 

®,® 

©-© 

©-© 


(D-© 

C 

.463 

0 


©-© 



©-© 

®-© 


— 

E 

.694 

A 

®-Q) 

©-© 

(D-0 

— 

— 

— 

— 

— 

F 

.925 

0- 

0-® 

©O 

Q-® 

— 

— 

— 

— 

— 


^■Flow configuration; 

@ Tunnel clean /S, 

@ Boundary- layer trip at nozzle entrance A 
(2) Rough side walls to within 40 inches of observing window 
Rough side walls to within 10 inches of observing window 
shock-induced separation 
© ©■•■ shock-induced separation |/^ 



A 

A 
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TABLE III.- RANGE OF OPTICAL VARIABLES 


Angle of 
incidence 
cp, deg 

No. of 
trans- 
missions, 
n/2 

Focal 

ratio 

f/no. 

Filter 

desig- 

nation 

Wratten 

no. 

Angular 

aperturesf* 

% 

(table I) 

Aerodynaml 

,c variables, table II 

Boundary- 

layer 

configu- 

ration 

Mach no. 

Density- 

ratio 

0 

1, 2 

8.6 

None 

1- 

-10 

(D-® 

V-VIII 

A 

0 

1, 2 

8.6 

None 



©-® 

III, V-VIII 

B 

0 

1, 2 

8.6 

None 



©-© 

II,III,V,VI, 

c 








VII 


0 

1, 2 

8.6 

None 



©-© 

II, III 

E 

0 

1, 2 

8.6 

None 



©-© 

I-III 

F 

■bij.5 

2 

8.6 

None 



© 

III,V 

B, C 


2 

8.6 

None 



© 

III,V 

B, C 

0 

1, 2 

12.3 

None 



©,© 

VIII 

B 

0 

1, 2 

17.2 

None 



©,© 

VIII 

B 

0 

OJ 

•N 

1 1 

^ 28.6 

None 



© 

III 

F 

0 

1, 2 

^ 2.9 

None 



© 

III 

F 

0 

2 

%.8 

None 



© 

VII 

C 

0 

2 

8.6 

45 



© 

III 

F 

0 

2 

8.6 

61 

\ 

! 

© 

III 

F 


^Series "b” apertures used for configuration only 
^Special plotting symbol for cp = 45 ^: III - C: 


V - C: O 


^Entrance pupil 
3/4 inch off 
^Entrance pupil 
1/2 inch off 
^Entrance pupil 
1/4 inch off 


covered by opaque mask with 
optical axis 

covered by opaque mask with 
optical axis 

covered by opaque mask with 
optical axis 


3 / 4 - inch-diameter hole 
1/2-inch-diameter hole 
1/ 4 - inch- diameter hole 


centered 

centered 

centered 
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TABLE IV.- TABULATION OF DATA 


E^^/Eq for aperture diameters, 
dxlO^ Inch 

^ .09 I ^*.86 I 4.20 I 3.32 j 1.80 f 1.28 I 1.01 I 0.75' 


(a) Boundary- layer configuration (^; cp ° 0°, D = 2.^ In. 


%/®o aperture diameters, 
liich 


8.09 I 4.86 I 4.20 I 3.^3^! 1!^^ 1.28 I 1.01 i 0, 


(b) Boxmdary-layer configuration (2 


0. UO 

.69 

.68 

.69 

.66 

.69 

.68 

.69 

.66 

.82 

.82 

.81 

.81 

.82 

.82 

.81 

.81 

1.43 

1.1^3 

1. U2 
l.kk 
1.^3 
1.42 

1.78 

1.78 

1.77 

1.75 

1.76 
1.76 

1.98 
2.00 
2.00 
2.04 
2.00 

1.99 
2.45 
2.50 
2.45 
2.45 


0.928 

.228 

.457 

.680 

.922 

.228 

.457 

.680 

.922 

.228 

.456 

.693 

.920 

.228 

.456 

.693 

.920 

.112 

.226 

.452 

.116 

.221 

.452 

.118 

.224 

.444 

.115 

.227 

.448 

.112 

.232 

.463 

.116 

.230 

.465 

.114 

.224 

.114 

.224 


458 

316 

583 

913 

1110 

316 

583 

913 

1110 

412 

824 

1160 

1530 

412 

824 

1160 

1530 

926 

1740 

3100 

928 

1730 

3150 

1380 

2270 

3960 

1140 

2260 

4230 

1370 

2700 

4760 

1410 

2570 

5020 

2150 

3790 

2280 

3950 


0.09 

.12 

.12 

.11 

.11 

.12 

.12 

.11 

.11 

.U 

.12 

.11 

.11 

.11 

.12 

.11 

.11 

.15 

.15 

.17 

.15 

.15 

.17 

.16 

.18 

.16 

.16 

.18 

.16 

.15 

.14 

.14 

.15 

.14 

.14 

.15 

.14 

.15 

.15 


1.00 

1.00 

1.00 

.999 

.996 

1.00 

.999 

.994 

.990 

1.00 

.999 

.996 

.989 

1.00 

.995 

.986 

.976 

1.00 

.998 

.988 

.999 

.994 

.976 

1.00 

.996 

.982 

.997 

.989 

.963 

.998 

.991 

.972 

.984 

.997 

.936 

.998 

.986 

.995 

.976 


1.00 

1.00 

.997 

.999 

.992 

.994 

.993 

.986 

.976 

1.00 

.999 

.990 

.981 

.996 

.985 

.971 

.941 

.997 

.994 

.969 

.999 

.987 

.951 

1.00 

.988 

.957 

.987 

.969 

.914 

.996 

.981 

.931 

.966 

.987 

.853 

.989 

.973 

.984 

.938 


1.00 

.995 

.994 

.999 

.991 

1.00 

1.00 

.990 

.980 

1.00 

1.00 

.986 

.968 

1.00 

.990 

.965 

.927 

.997 

.998 

.962 

1.00 

.989 

.952 

1.00 

.986 

.949 

.989 

.968 

.892 

.992 

.974 

.908 

.961 

.987 

.816 

.989 

.970 

.990 

.925 


1.00 

1.00 

1.00 

1.00 

.981 

.995 

.990 

.982 

.961 

.996 

.994 

.982 

.950 

.996 

.975 

.933 

.888 

.998 

.987 

.965 

.995 

.972 

.903 

.995 

.976 

.933 

.976 

.952 

.833 

.971 

.962 

.860 

.929 

.977 

.720 

.989 

.955 

.979 

.886 


1.00 

1.00 

.980 

.989 

.953 

.990 

.964 

.919 

.875 

.995 

.981 

.861 

.980 

.921 

.818 

.685 

.985 

.960 

.865 

.985 

.927 

.733 

.980 

.9'+9 

.780 

.957 

.869 

.595 

.968 

.886 

.656 

.801 

.937 

.441 

.963 

.847 

.918. 

.706 


0.992 

.991 

.981 

.954 

.902 

.962 

.928 

.837 

.769 

.984 

.947 

.870 

.761 

.960 

.862 

.710 

.540 

.971* 

.936 

.753 

.965 

.853 

.557 

.964 

.895 

.639 

.917 

.764 

.417 

.93** 

.803 

.498 

.667 

.903 

.286 

.932 

.739 

.864 

.564 


0.943 

.984 

.961 

.906 

.836 

.952 

.897 

.755 

.673 

.976 

.907 

.807 

.672 

.931 

.786 

.599 

.421 

.946 

.884 

.625 

.910 

.762 

.418 

.941 

.837 

.520 

.872 

.661 

.306 

.909 

.712 

.378 

.526 

.847 

•195 

.863 

.598 

.786 

.441 


0.920 

.968 

.933 

.874 

.763 

.945 

.864 

.686 

.592 

.973 

.855 

.731 

.578 

.925 

.709 

.503 

.330 

.931 

.838 

.535 

.890 

.671 

.336 

.935 

.762 

.416 

.790 

.569 

.218 

.877 

.613 

.267 

.460 

.784 

.112 

.816 

.515 

.714 


0.41 

.65 

.65 

.65 

.64 

.68 

.68 

.68 

.67 

.80 

.80 

.79 

.79 

.82 

.81 

.80 

.81 

1.43 

1.43 

1.42 

1.43 

1.42 

1.42 

1.76 

1.75 

1.76 
1.76 
1.76 
1.75 
2.01 
2.01 

1.96 
2.00 
1.99 

1.97 

2.44 
2.50 
2.44 


.335 2.51 


0.936 

5 780 

(1) 

.236 

; 560 


.465 

1150 

— 

.701^ 

1510 

— 

.936 

2080 

— 

.232 

630 

— 

.457 

1210 



.686 

1810 



.921 

2180 



.233 

810 



.465 

1600 

— 

.703 

2220 



.936 

2940 



.231 

750 

— 

.463 

1470 

— 

.696 

1970 



.922 

2780 

— 

.112 

1120 

.33 

.225 

2190 

.25 

.1*53 

3820 

.22 

.115 

1370 

.33 

.235 

2330 

.25 

.1*53 

3910 

.22 

.111 

1420 

.23 

.228 

2560 

.22 

.461 

531*0 

.23 

.113 

1390 

.23 

.223 

2610 

.22 

.461 

5340 

.23 

.109 

1660 

.20 

.225 

2970 

.17 

.474 

51*70 

.18 

.110 

1640 

.20 

.228 

3060 

.17 

.474 

5790 

.18 

.115 

2360 

.19 

.225 

4020 

.18 

.114 

2290 

.19 

.221 

4100 

.18 


2 

1.00 

1.00 

0.996 

1 1.00 

1.00 

0.975 

1 

1.00 

1.00 

.99c 

» 1.00 

1.00 

1.00 

1 

1.00 

1.00 

1.00 

1.00 

1.00 

.982 

1 

.998 

1.00 

.992 

.999 

.97^ 

^ .941 

1 

.996 

.996 

.991 

.988 

.951 

.893 

2 

1.00 

1.00 

1.00 

1.00 

1.00 

.994 

2 

.998 

1.00 

.997 

1.00 

.969 

.929 

2 

.994 

.997 

.994 

.992 

.927 

.845 

2 

.983 

.982 

.973 

.964 

.847 

.711 

1 

1.00 

.998 

.994 

.995 

1.00 

.995 

1 

.999 

1.00 

1.00 

1.00 

.970 

.939 

1 

.995 

.985 

.992 

.983 

.947 

.871 

1 

.987 

.980 

.975 

.964 

.864 

.752 

2 

.999 

.987 

.999 

.989 

.992 

.969 

2 

.991* 

.988 

1.00 

.989 

.958 

.872 

2 

.984 

.971* 

.977 

.962 

.868 

.705 

2 

.973 

.957 

.958 

.915 

.754 

.559 

1 

1.00 

.995 

.995 

.995 

.984 

.960 

1 

.998 

.987 

.981 

.974 

.952 

.901 

1 

.989 

.978 

.967 

.965 

.880 

.767 

2 

.985 

.973 

.968 

.965 

.960 

.940 

2 

.993 

.977 

.976 

.960 

.912 

.860 

2 

.976 

.938 

.927 

.888 

.713 

.551* 

1 

1.00 

1.00 

1.00 

.993 

.970 

.950 

1 

.991 

.975 

.973 

.975 

.917 

.808 

1 

.968 

.928 

.907 

.881 

.682 

.481 

2 

.996 

.992 

.988 

.971* 

.948 

.892 

2 

.985 

.957 

.963 

.945 

.828 

.708 

2 

.945 

.858 

.834 

.756 

.465 

.306 

1 

.998 

.993 

.988 

.988 

.971 

.936 

1 

.987 

.977 

.971 

.959 

.891 

.762 

1 

.958 

.911 

.885 

.816 

.573 

.396 

2 

.995 

.990 

.977 

.978 

.937 

.886 

2 

.980 

.958 

.947 

.909 

.753 

.611 

2 

.931 

.827 

.770 

.677 

.376 

.246 
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212 


eirithmlc velocity decrement law: 


I were given to good approximation by the log- 


/ \i/m 


where 6 in this case is the half -width of the wind tunnel (6 in.). The value of m for all the profiles fell in the range 

5.6 < m < 6.2 

The lower value is associated with the lowest density ratio (p^g = 0.22). 
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(a) Tmnel Mach n™ber = 0; missile Mach number = h 



(b) Tunnel Mach number = 2; missile Mach number = 5 

Figure 1.- Shadow photographs of missile models in the Ames supersonic 
free-flight wind tunnel showing the effect of turbulent layers on a 
collimated light beam; spark duration O .5 microsecond. 
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Figure 2.- Experimental arrangement for viewing a point light source at infinity through two or 
four thickened turbulent wall boundary layers in the Ames 1- by supersonic wind tunnel 

no. 1. 


NACA RM A56B21 


MCA RM A56B21 


41 




Figure 3*" Variation of the fraction of radiant flux received through 

four turbulent boundary layers with the fraction received through two. 
Boundary- layer configuration (See table I for symbol identification.) 
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(a) Effect of changing houndary-layer thickness and shape at constant free-stream Mach ntimher and 

free- stream density. 


Figure 4.- Variation of the fraction of radiant flux received through two turbulent boundary layers 

with the angular aperture of the telephotometer. 
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(b) Correlation of data for various Mach numbers and densities employing parameters predicted by 

scattering theory. (See table II for symbol identification.) 


Figure 4.- Concluded. 
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15 sec 



(a) 



Exposure time 
45 sec 
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Figure 5 .- Comparison of relative intensity distributions with photo- 
micrographs (X320) of the image at the focal plane of the receiving 
telescope for various wind-tunnel flow conditions. Light transmis- 
sion through two boundary layers; boundary-layer configuration(^. 
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Figure 6 .- Effect of various meihods for "thickening ■turhuleni boundary layers on "the shape of "the 

Mach number profiles. 
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Figure 7 *- Attenuation of radiant flux through one, two, and four turbulent boundary layers as a 
function of the combined mean flow and optical parameter, ^/e^cos cp. (See tables II and III 
for symbol identification.) 
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